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Abstract 



We present an exact, axially symmetric, static, vacuum solution for f(R) gravity in Weyl's canonical coordinates. We obtain a 
general explicit expression for the dependence of df(R)/dR upon the r and z coordinates and then the corresponding explicit form 
of f(R), which must be consistent with the field equations. We analyze in detail the modified Schwarzschild solution in prolate 
spheroidal coordinates. Finally, we study the curvature invariants and show that, in the case of f(R) + R, this solution corresponds 
to a naked singularity. 
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Q 1. Introduction 

In recent years, f(R) theories of gravity have gained much at- 
r ^ 1, tention as promising candidates to overcome the issues posted 
Q-iby dark energy in the standard cosmological model (c.f. refer- 
^ -ence [1, 2] for a recent review). There has been a stimulating 
(^debate in their study, leading us to a number of interesting re- 
'— '-suits, particularly in the context of exact solutions. 

| Due to the highly non-linear nature of the /(i?)-field equa- 
-tions, finding exact solutions is indeed a difficult task. Fol- 
ly") lowing astrophysical motivations, solutions with spherical sym- 
OO metry have been the most widely studied [3]. However, there 
(N has also been a growing interest in finding exact solutions with 
^1 cylindrical symmetry [4]. To the best of the authors' knowledge 
^ — 1 (except for [5, 6]), there is no fully integrated, explicit and exact 
axially symmetric solution of /(/?)-gravity. 

In this paper, we consider static vacuum solutions of f(R) 
• • theories in Weyl coordinates. In particular, we obtain the ex- 
# ^ plicit dependence of df(R)/dR upon the coordinates p and z. 
-This in turn, allows us to get the corresponding explicit form of 
$— i f(R)- Finally, we analyze in detail the solutions of the modified 
^ field equations corresponding to the Schwarzschild solution in 
cylindrical coordinates. 

2. /(/?) Field equations for a static axially symmetric 
space-time 

The f(R) action is given by 



(N 



s = I[ihG f{R)+ ^)^ d4x > 



(1) 



where G is the gravitational constant, R is the curvature scalar 
and £ m is the matter Lagrangian. The field equation resulting 
from this action are 



>ab 



-Rg ah = %nG{f[ \ + T ™), 



where the 'gravitational' stress-energy tensor is 



8jtGT s , = T\ , 

ab ab ' 



(2) 



(3) 



and 



ab F R 



^ (f(R) - RF R ) + V c V d F R [5<6 b d - g ab g cd ) 



Here F R = df(R)/dR and = T™/F R , where is the stress- 
energy tensor obtained from the matter Lagrangian £ m in the 
action (1) . 

Equivalently, we can write (2) in the form 

F R R ab - \f{R)g ah - V a V 6 F fl + gab UF R = 8nGT™. (4) 

Taking the trace of this expression we obtain the relation be- 
tween f(R) and its derivative F R 



F R R - 2f(R) + 3DF R = %nGT u 



(5) 



We are interested in the static axially symmetric solutions of 
(4). To this end, let us consider the Weyl- Lewis-Papapetrou 
metric in cylindrical coordinates is [7] 



ds 2 = -e^dt 2 + e- 2 ^\p-dur + e 2A (dp 2 + dz 2 )], 



(6) 
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where cp and A are continuous functions of p and z. Using the 
trace equation (5), the modified Einstein field equations (4) be- 
come 



F R R ab - V fl V feJ F R - 87rGr™ = g ab B 



(7) 
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where B = j(FrR — DF r - %nGT m ). The non-zero components 
of the Ricci tensor are 



3. /(/?) vacuum solutions for a static axially symmetric 
space-time 



^00 


= e 4 *- u V 2 cf>, 


(8a) 


Rn 


= pV 2i V 2 0, 


(8b) 


R22 


= -V 2 A + V 2 + -A, p - 2(0, p ) 2 , 
P 


(8c) 


R33 


= -V 2 A + V 2 - 2<f>l , 


(8d) 


R23 


= ~ A ,z ~ 2 <f>,p<P,z, 

P 


(8e) 



while a straightforward computation of the curvature scalar 
yields 



For simplicity, we restrict ourselves to the vacuum case, i.e. 
we make 



0. 



b 2 : ) + z£r(F Rj)p - F R , ZZ ), 

Zf R 



A tZ = 2p4> :P <f) 7Z + 



pF R , pz 



pF R , z (F R , pp - F RtZZ ) + P F R V 2 (F R , Z ) 

+ F R , zp (F R -2pF R , p ) = 0, 

pp-< < v, p ■ r, z / 



R 



f(R) = IfjjR + 3 -UF R . 



(12a) 
(12b) 

(12c) 

(12d) 
(12e) 

(12f) 



R = 2e 2 ^ 2 ' 1 (-V 2 i + V 2 + —A 



i.2 ±2 



), (9) 



where V 2 is the usual Laplace operator in cylindrical coordi- 
nates. 

From (7) and (8) we obtain the following system of equa- 
tions: 



— \F R R m - inGT^] = B, 

goo 

— [F R R n -8nGT™]=B, 
gn 

— [F R R2~> - ZnGT^,] = B, 

g22 

— \f r R 3 s - %nGT™A = B, 

#33 L J 

FrRt3 - SnGT^ - F R 23 = 0. 



(10a) 

(10b) 

(10c) 

(10d) 
(10e) 



This allows us to write down the independent field equations, 



i.e. 



V 2 0= _^ e 2(A- W[r mO_ r ml ]) 

F R 

A_ n = p(4, 2 n - fa + ^E e ^\ T f - Tf] 



> ^ P(<P, P ■ 

P 
2F f 



F R 

[F R ,22 - F R ^], 



] O A, A, ^ 87rG P -rm ^ P F R,23 
\z = ^P<P,p<f>,z + —=T- T 23 + S, — ■ 



(11a) 

(lib) 
(11c) 



One can see that equation (12d) is the integrability condition 
for A. Note that for an arbitrary f(R), the system (12) may be- 
come inconsistent. Therefore, we look for the class of functions 
f(R) compatible with (12). It is an easy exercise to prove that 

DF R = e 2( *- A \F R , pp + F RiZZ ). (13) 

Thus, substituting (13) in (12f) and using (12e) we obtain 
3W(p) 



f(R) = -^F R R 



1 - 



2F R 



where 
W(p) 



Fr, pp + F R 



(14) 



(15) 



In order to obtain some analytical solutions to the system 
(12), we need to make some further simplifying assumptions. 
First, suppose that it is possible to write 



F R (p,z) = U(p)V(z). 
Then, substituting back into (12d) we have 



(16) 



P - l u- 2 



2 dU I dU jj\ 2 jj d2u 
dp \ dp j dp 2 



dVY 



d 3 V d 2 VdV 
dz 3 dz 2 dz 



(17) 



Equating each side to a separation constant, I 2 , we obtain the 
third order pair of ordinary differential equations 



p- x u- 2 



2 dU I dU \ 2 v d 2 U 
dp \ dp j dp 2 



*2 



U y d 3 V d 2 VdV 
dz 3 dz 2 dz 



= l 2 



I 2 . 



(18a) 
(18b) 



To find a general solution to the above equations is indeed 
a difficult task. Nevertheless, in the following sections we will 
discuss some particular solutions to (11a). 



Let us re- write (18a) as 



dM(p) M(p) 
dp p 



'2' 



(19) 



where M(p) = U l dU/dp. One can solve this immediately to 
obtain 



i dU n l 2 p 
M(p) = U~ l — = -~^-, 
dp p A 

which has the solution 

U(p) = cp"e-' lpl/ \ 

where c and n are integration constants. 
Now, one can easily show that 

V(z) = e bz , 



(20) 



(21) 



(22) 



where b is an arbitrary constant, is solution of (18b) if both, b 
and /, satisfy the condition 



bl z = 0. 



(23) 



Thus, we have that some possible solutions for Fr are [c.f. 
equation (16)] 



(i) b — and / = 0. In this case 

Fr = cp". 

(ii) b + and I — 0. In this case 

F R = cp"e b: . 

(iii) b — and / + 0. In this case 



Fr = cp' 



(24) 
(25) 
(26) 



Consequently, by substituting (24) [or (25)] in (14) and using 
(12) we obtain that f(R) must satisfy the consistency condition 



df 

f(R) = 2R—, 
JK ' dR 

whose solution is simply 



f(R) = kR ll \ 



(27) 



(28) 



where k is an arbitrary constant. 

Substituting (25) in both, equations (12b) and (12c), we have 



A p = p(cf> 2 p - 2 ) + — [n(n -l)-bYl (29a) 



1 



^2„2 n 



2p 



A- 



> , + bn, 



(29b) 



respectively. Whereas, by taking Fr = cp"e hp '^ we obtain 
from (14) 



f(R) = ±F R R [1 - 3L(p)] , 



with 



Up) 



(fp 2 - An) 2 - 4(l 2 p 2 + Ari) 



(3Pp 2 + An - A)(l 2 p 2 - An) ' 
Finally, substituting (26) in (12b) and (12c), we have 

+ -L[(/ 2 p 2 -4n) 2 -4(/ 2 p 2 +4n)], 
32p 



(30) 



(31) 



(32a) 
(32b) 



respectively. 

As we can see from equations (29) and (32), the function A 
can be calculated by means of a line integral. Although V 2 <p = 
is a linear differential equation, the equations for A manifest the 
non-linearity of the "modified" Einstein field equations. 

The usual Einstein vacuum equations for the static axisym- 
metric spacetime we have 



V 2 



0. 



A fi = P(<P, P - <PJ> 
A^ z — 2p(p p (p a . 



(33a) 
(33b) 
(33c) 



The Laplace equation may be solved by using various coordi- 
nates in the Euclidean 3-space and then the function A can be 
calculated (c.f. chapter 20 in [7]). We observe the following 
points: 

1. If Fr = cp n e bz , then we can obtain a vacuum static axially 
symmetric solution (cp, A) of the vacuum 'modified' Ein- 
stein field equations from the vacuum Einstein field equa- 
tions 0, A) using the transformation 



A = X + In [/tp" ( "- 1)/2 ] p 2 + bnz. 



2. Similarly, if Fr = cp"e '" p "^ 8 we make 



A = A 



(34a) 
(34b) 



(35a) 



+ 2(2,1 + 1)lV + ln K 6 ""^]} • < 35b > 

In both cases, k is an appropriate constant. 

Thus, one can say that equations (34) and (35) represent a 
Weyl class of solutions in /(/?)-gravity. Moreover, one can ex- 
pect that some properties of the curvature of the seed solution 
will be inherited to the modified ones. 



4. A particular solution for a vacuum static axially sym- 
metric space-time 

Here we present an application of the results obtained in the 
previous section. First, we assume a given metric potential, <p 
say [c.f. equation (6)], and then we find the other one by means 
of the two equations Fr = cp"e hz , and Fr = cp"e~'~ p ~l & . Let us 
work in prolate spheroidal coordinates (x,y), with x € [l,°o) 
and y € [—1,1]. These are related to the cylindrical coordinates 
(p, z) through the relations 



p 2 — m 2 (x 2 - 1)(1 — y 2 ) and z — mxy. 



(36) 



The line element (6) becomes (c.f. equation (4.5.18) in [8]) 

dx 2 dy 2 



j 2 2d> A J. , 2 2(A-(h) I 2 2\ 

as = - e v at + m e y [x - y J 



x 2 -l l-y 2 



(37) 



The Einstein vacuum equations in f(R) gravity for a static 
axially symmetric space-time can be cast into the form 



where 

<P-- 



ln 





A,x + Pp,x + &Z,x, 



and A 



x+l 



In 



x 2 — y 2 



(38a) 
(38b) 
(38c) 



(39) 



are the metric potentials corresponding to the usual 
Schwarzschild solution in standard general relativity. Note 
that we can relate the prolate coordinates (x,y) to the usual 
Schwarszchild coordinates (r, 6) through the relations (c.f. 
equation (4.5.19) in [8]) 

r = m(x +1) and 9 = arccos(y). (40) 

This in turn allow us to relate them to (36) by 

p 2 = [(r - to) 2 - m 2 ] sin 2 (#) and z = (r — m) cos(fl). (41) 

It becomes clear that p — corresponds to x = 1 and r = 2m. 
Thus, in the forthcoming discussion, the reader should be aware 
that the domain of the prolate coordinates is defined from the 
horizon up to infinity. Notice as well that 



P - ^t( f R,pp 



F RtZZ ) and Q. 



PF R , PZ 



(42) 



Thus, in the same way as the last case, we will obtain a explicit 
form of A by considering the different values of F R . 



4.1. F R = cp"e bz 
Here we have 



n(n — l)x 

(x 2 - l)(x 2 -y 2 ) *" 2(x 2 - 1) 
b 2 m 2 



y 



-x(\ — y ) + bmny 
n{n — l)y 



x 2 -y 2 2(1 -y 2 ) 
b 2 m 2 y 2 

— - — {x — 1) + bmnx, 



(43a) 



(43b) 



whose solution is 
A 



^^-^Inf^-lXl-y 2 )] 



Q 



b 2 m 2 Q 



2 2 

xr y 



+ bmnxy, 



■y 



(44) 



(45) 



Substituting (39) and (44) into the metric (37), we observe 
that there are some apparent singularities whenever x and y take 
the values ± 1 . Note that x — - 1 is not part of the domain of the 
prolate coordinate system and that the singularities in y are the 
trivial pole singularities when 6 = and 6 — n. The singularity 
in x — 1, however, deserves some deeper analysis. To this end, 



we compute the two main curvature invariants, i.e. R = R° and 

A - KabcdK 

Let us assume that n > 1 and consider the two possibilities, 
b = and b + 0, separately. 



1 . b — 0. In this case, the Ricci curvature scalar is 



R = - 



n(n-l)(x 2 -y 2 )(x 



2 2 2 



l + y 2 ) 



-l/2n(n-l) 



m 2 (x+ 1) 2 (-1 +y 2 )(x 2 - 1) 



whereas the Kretchman invariant takes the form 

f(x,y,n) 



K = 



2 (x- \) 2 m A (x+ 1) 6 (1 -y 2 ) 2 ' 



(46) 



(47) 



where /(x, y; n) is a lesser illuminating expression. One 
can easily show that in the large x limit, both of the above 
expressions converge to zero. Moreover, inspecting the 
« = 1 solution in the 'equatorial' plane, y = 0, one obtains 



and 



K„ = \ — 



Rn=l = 



44 44m 2 



m 4 (x+l) 6 r 1 



.6 ' 



(48) 



(49) 



in agreement with the Schwarszchild solution. 

2. b + 0. This case is richer in content. The curvature invari- 
ants can be written as 



R 



and 



K=-- 



h\(x,y; n) 
m 2 (x +1) 

h 2 (x,y; n) 



2{x- l) 2 m 4 (x+ 1) 6 (1 -y 2 ) 



2 ' 



(50) 



(51) 



Just as before, the functions /i,(x, y; n) are long polynomial 
expressions from which little can be said. Interestingly, a 
similar analysis of the « = 1 soution in the y = plane and 
taking b — 1 yields 



Rn=\ — —■ 



and 



K„ = \ — 



J l/2m 2 x 2 x 2 

(x+lf 
Oix 1 ) 



(x+lf(x- \)m A 



(52) 



(53) 



Here we observe that the regular behaviour of the n = 1 
solution at the 'horizon' of the b — case is lost. What 
we see here is a true singularity at x = 1 with no horizon 
dressing it. 

The behaviour in the vicinity of x — 1 is shown in figure 1 
for the soutions with n = 1, 2, 3. 



4 



n=l n=2 n=3 





n=l n=2 n=3 




Figure 1 : Kretchman invariant (top plots) and Ricci curvature scalar (bottom plots) as functions of x for Fr = cp" e bz with n = 1 , 2 and n = 3 for the values b = 
(left plots) and b = 1 (right plots). 
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n=l n= 2 n= 3 



n = 1 n=2 — ■ — n = 3 





10-f 



R 



- 10-1 



n=l n = 2 n=3 



R 



1.5 



2.5 




-30-1 



2 2.5 



Figure 2: Kretchman invariant (top plots) and Ricci curvature scalar (bottom plots) as functions of x for Fr = cp"e '"P"/ 8 with n = 1,2 and n = 3 for the values 
/ = (left plots) and I = 1 (right plots). 
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4.2. F R = cp"e-' 2 P 2/8 
Now, in this case 



A . = A x + z P, 

• 32(x 2 -l) 



Ay = Ay- 



y 



-p. 



32(1 -y 2 ) 

P = [16«(«-l) + Z 4 p 4 -4/ 2 (2« + l)p 2 ], 
and the solution is written as 

A = A + In [(x 2 - l)(l-y 2 )] 

+ ^[/ 2 m 2 (e + 2) + 8(2«+l)]e. 

Just as in the previous case, we analyse the curvature invari- 
ants to look for singularities. Again, we split our study into two 
cases 

1 . I — 0. The curvature invariants are 

h 3 (x,y;n) 



R 



and 



K 



(x + l)m 2 
h\{x, y; n) 



2 (jc- \fm A {x+\) b (-1 +y 2 Y 



static and axially symmetric vacuum solutions of f(R) + R cor- 
respond to naked singularities, as can be seen in the right hand 
column in figures 1 and 2. In particular, one observes that the 
(54a) singularity structure of the case Fr = cp"e hz is very sensitive 
to the value of b as can be seen in figure 1 where the solution 
(54b) ceases to be regular at the horizon and becomes a naked singu- 
larity, even in te n — 1 case. 

Finally, it is worth noting that the potentials (44) and (55) 
were found by integrating the corresponding system of differ- 
ential equations [equations (38)]. However, these solutions can 
be obtained directly by using the transformations (34) and (35). 
This result allows us to generate axially symmetric solutions 
for f(R) from known seeds of the Weyl class of solutions for 

(55) 

the Einstein field equations. 

This work is dedicated with great pleasure to Biky 
(M. V.R.H.) on the occasion of her 23rd birthday. 
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Same as before, the equatorial plane n — 1 solution re- 
duces to Schwarszchild 



Rn=\ = 



and 



K n= \ — 32 ■ 



1 



(58) 



(59) 



(x+ l) 6 m 4 ' 

2. I + 0. This case shares the same singular structure as the 
I = case. As can be seen in figure 2 for n = 1,2,3. 
Note that for these soutions, n — 1 is always regular at the 
horizon. However, taking « > 1 always produces a naked 
singularity at x — I. 

Finally, let us note that, in every case, the curvature invariants 
converge to zero in large x limit. 
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5. Closing remarks 

The issue of static and axially symmetric solutions in f(R)- 
gravity is a timely topic in the context of the exact solutions. 
In this paper, we have presented an axially symmetric static va- 
cuum solution in Weyl coordinates for f(R) gravity. In particu- 
lar, from the integrability condition of one of the metric poten- 
tials of the Weyl-Lewis-Papapetrou line element and using the 
method of separation of variables we have obtained a general 
explicit expression for the dependence of df(R)/dR on the p 
and z coordinates and, therefore, the corresponding general ex- 
plicit form of f(R). Working in prolate spheroidal coordinates, 
we have analysed in detail the 'Schwarzschild' solution to the 
modified field equations. We have shown that these particular 
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